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Using ab initio methods, we have studied the electron-phonon interaction in compressed, body- 
centered cubic (bcc) yttrium, which is predicted to be stable at 280 GPa [Melsen et al, Phys. Rev. 
B 48, 15574 (1993)]. We find that compressed, bcc yttrium has a large electron-phonon coupling 
with A = 1.8, leading to a superconducting transition temperature T c m 30 K or above. Our results 
indicate that the large electron-phonon coupling is due to the lattice hardening. 



Recent experiments on pressure- induced increase in su- 
perconducting transition temperature T c of yttrium [l| 
seem to suggest an almost continuous increase in T c with 
pressure. With increase in pressure from ambient to 
115 GPa, the structural changes @, 0| in yttrium from 
hcp^Sm-type^dhcp^trigonal are accompanied by in- 
crease in T c from 6 mK Q to 20 K In contrast, T c 
increases with pressure only in the fee phase of lithium as 
it changes from bcc— >fcc^hRl^cI16 with pressure going 
up to 70 GPa Thus, the pressure-induced changes in 
T c seem to be intimately connected with the structural 
phase transitions in these metals [fj]. 

It has been shown theoretically that yttrium would 
stabilize in a bcc phase above 280 GPa Q. Taking a cue 
from the experiments on yttrium [H, it is not surprising 
to expect bcc yttrium to superconduct well above 20 K 
at pressures above 280 GPa. However, it is also possible 
that in this pressure range yttrium may start to behave 
like lithium H, @, ©, El | . 

To further investigate the changes in the supercon- 
ducting properties of bcc yttrium, we have used density- 
functional-based methods to calculate its electronic 
structure, lattice-dynamical response and the electron- 
phonon interaction. Subsequently, we have solved the 
isotropic Eliashberg gap equation to obtain the super- 
conducting transition temperature T c . 

Based on our calculations, we find that compressed, 
bcc yttrium has a large electron-phonon coupling with 
A = 1.8, leading to a superconducting transition temper- 
ature T c m 36 K for //* = 0.25. Such a large coupling 
is a direct consequence of lattice hardening and its sub- 
sequent coupling to Fermi electrons. We also find the 
electrons at Fermi energy to have predominantly ti g sym- 
metry. Before we describe our results, we briefly outline 
the computational details of our calculations. Further 
details can be found in our previous work [12] as well as 
the references therein. The notations used here are that 
of our previous work [H]- 

We have calculated the electronic structure, the 
phonon density of states, F(u>), and the Eliashberg func- 
tion, a 2 F(u>) of bcc yttrium using full-potential, lin- 
ear muffin-tin orbital (LMTO) method [H, 0]. Subse- 
quently, we have numerically solved the isotropic Eliash- 
berg gap equation [H, 0, E3] for a range of /i* to obtain 



the corresponding T c . 

The lattice constant a = 5.15 a.u. used in our calcu- 
lations for bcc yttrium approximately corresponds to a 
compressed volume equal to 0.31wo, where vq is the exper- 
imentally determined equilibrium volume of hep yttrium 

The charge self-consistent, full-potential, LMTO calcu- 
lations for bcc yttrium was carried out with 2k— energy 
panel, 897 k-points in the irreducible wedge of the Bril- 
louin zone, and the muffin-tin radius for yttrium was 
taken to be 2.22 atomic units. The 4s state of yttrium 
was treated as a semicore state, while the 4p state was 
treated as a valence state. The basis set used consisted of 
s, p, and d orbitals at the yttrium site, and the potential 
and the wave function were expanded up to l ma x = 6. 

The Fermi surface, phonon density of states, and the 
Eliashberg function of bcc yttrium were calculated us- 
ing the linear response code based on the full-potential, 
LMTO method. The Fermi surface was constructed using 
XCrySDen program [18] with eigenvalues calculated on a 
48 x 48 x 48 grid in the reciprocal space. The dynamical 
matrices and the electron-phonon matrix elements were 
calculated on a 8 x 8 x 8 grid resulting in 29 irreducible 
q-points. The Brillouin zone integrations during linear 
response calculations were carried out using a 16 x 16 x 16 
grid of k-points. The Fermi surface sampling for the eval- 
uation of the electron-phonon matrix elements was done 
using a 48 x 48 x 48 grid. Another set of calculations 
of self-consistent charge density, F(u>) and a 2 F(uj) was 
carried out at a lattice constant o = 5.1987 a.u. using a 
6x6x6 grid with 16 irreducible q-points. In this 
12 x 12 x 12 grid of k-points was used for linear-response 
calculations, in conjunction with a 36 x 36 x 36 grid for 
Fermi surface sampling. The muffin-tin radius for yt- 
trium was taken to be 2.251 atomic units. For visualizing 
the charge density within a 5 mRy energy window to the 
Fermi energy we used the atomic-sphere-approximation 
variant of the LMTO code [H]. 

In Fig. 1, we show the total and the /-resolved densities 
of states (DOS) of bcc yttrium. For comparison, we have 
also included the DOS for the hep yttrium [12j. From 
Fig. 1, we see that the pressure-induced pushing out of s- 
electrons is followed by creation of new d-states between 
—0.55 Ry and —0.25 Ry to accommodate the s-electrons. 
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Figure 1: The s, p, d and the total densities of states of 
bcc (solid line) and hep (dashed line) yttrium. The vertical 
dashed line indicates the Fermi energy. 
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Figure 2: The symmetry-resolved d densities of states in bcc 
(top panel) and hep (bottom panel) yttrium. The densities of 
states in the top panel correspond to t^ g (xy, yz and zx; solid 
line) and e g (x 2 — y 2 and 3z 2 — 1; dashed line) symmetries, while 
in the bottom panel they correspond to yz (zx; solid line), xy 
(x 2 — y 2 ; dashed line) and 3z 2 — 1 (dotted line) symmetries. 
The vertical dashed line indicates the Fermi energy. 



The flattening of the rf-band reduces the DOS at the 
Fermi energy to 10.03 st/Ry from its hep value of 25.89 
st/Ry, with most of the contributions coming from the 
d-states. 

It is interesting to consider the symmetry-resolved d- 
DOS, as shown in Fig. 2. Some of the differences in the 
DOS between the ambient pressure hep phase and the 
high-pressure bcc phase arise due to change in symmetry 
from hexagonal to cubic. In hep phase, the three inde- 
pendent rf-DOS corresponding to yz (zx), xy (x 2 — y 2 ) 
and 3z 2 — 1 symmetries have similar bandwidths but some 
differences in the DOS as a function of energy. However, 
in bcc phase the bottom of the doubly degenerate e g 
(x 2 — y 2 and 3z 2 — 1) symmetry DOS is almost 0.2 Ry 
below the bottom of the triply degenerate ti g (xy, yz and 




Figure 3: The band-structure of bcc yttrium (solid line) along 
high symmetry directions in the corresponding Brillouin zone. 
The horizontal line, passing through the energy zero, indicates 
the Fermi energy. 

zx) symmetry DOS, indicating transfer of pushed-out s- 
electrons to e g symmetry orbitals. A symmetry-resolved 
breakup of DOS at the Fermi energy is 0.40, 1.14, 5.87 
and 2.62 st/Ry for s, p, t2 g and e g symmetries, respec- 
tively. 

As a prelude to understanding the Fermi surface topol- 
ogy and the nature of electrons contributing to the 
electron-phonon interaction, we show in Fig. 3, the calcu- 
lated band-structure of bcc yttrium along high symmetry 
directions in the corresponding Brillouin zone. The band 
corresponding to 4p state is much lower in energy, and 
hence it is not shown in the figure. Based on Fig. 3 and 
an analysis of the eigenvectors (not shown in the figure) , 
the so-called fat bands, we find that the bands crossing 
the Fermi energy along T-N, N— P and H— P symmetry 
directions have substantial ti g symmetry states, indicat- 
ing an important role for these electrons as far as the 
electron-phonon interaction was concerned. 

Since Fermi surface plays an important role in under- 
standing many of the superconducting properties of ma- 
terials, we show the calculated Fermi surface of bcc yt- 
trium in Fig. 4. The Fermi surface due to band number 
1 surrounds T point and it has a bucket-like structure 
around P symmetry point in the reciprocal space. The 
second band has an umbrella-like structure with four-fold 
symmetry around N point, in addition to having surfaces 
around T and P points. It turns out that some of the flat 
portions in all these surfaces lead to Fermi surface nest- 
ing especially along N— P and T— P symmetry directions, 
the indications of which can be seen in Figs. 4 (c) and 
(d), showing the Fermi surface cut in a plane parallel to 
(110) plane. A more detailed analysis is needed to clearly 
identify all the nesting vectors. 

In Fig. 5 we show the charge density within a 5 
mRy energy window around the Fermi energy in bcc yt- 
trium, calculated usin g th e tight-binding-LMTO method 
as mentioned earlier [13]. The e g symmetry electrons 
seem to be localized around the yttrium sites while ti a 
electrons are more spread out and hybridized with the 
nearby yttrium atoms. Having described the electronic 
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Figure 4: The Fermi surface of bcc yttrium in reciprocal space 
unit cell corresponding to (a) band number 1, and (b) band 
number 2. The Fermi surface cut on a plane parallel to (110) 
plane for (c) band number 1, and (d) band number 2. The 
center of the cube corresponds to the F point. 
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Figure 5: The isosurfaces and contours of charge density 
within a 5 mRy energy window around the Fermi energy in 
bcc yttrium. The values of the isosurfaces are (a) 3.5 and (b) 
4.5 in units of 10 -2 [e/a.u. 3 ], which is also the unit used in (c) 
and (d). 



structure of bcc yttrium, we turn to its lattice-dynamical 
response and how the phonons couple to electrons close 
to Fermi energy. 

Our calculated phonon density of states F(u>) and the 
Eliashberg function a 2 F(uj) of bcc yttrium are shown in 
Fig. 6, where we also show the corresponding values for 
hep yttrium taken from our earlier work In the in- 
sets of Fig. 6, we show the variations of u> and A q „ of bcc 
yttrium along T-N direction. A comparison of F(u>) of 
bcc yttrium with that of hep yttrium clearly shows the 
hardening of the bcc lattice, which has become almost 
three times as hard. Some of the significant additions 
to F(uj) of bcc yttrium take place around the peak at 
59 meV, the peak itself arising out of the displacement 
of yttrium atoms with components in all the three direc- 
tions. We also find that the phonon mode (with displace- 
ment having components along +y and — z directions) 
in bcc yttrium along F-N direction softens to 9 meV at 
q = (0.0, 0.375, 0.375) in units of 2n/a, and it has a value 
of 12.6 meV at the N symmetry point. Phonon softening 
in fee yttrium is discussed by Yin et al. ficj ]. 

The pressure-induced increase in a 2 F(uj) in yttrium, 
accompanied by structural phase transformations, seems 
to arise from the enhancement of contributions from ex- 
isting modes and coupling strongly to the newly created 
modes. For example, the peaks around 23 and 27 meV 
are created in the dhep phase and they become much 
stronger in the bcc phase, while the peak at 59 meV 
is created in the bcc phase. As a result, the average 



electron-phonon coupling constant A increases from 0.46 
in hep yttrium to 1.8 in bcc yttrium. 

Somewhat surprisingly, we find that the largest contri- 
bution of A qi/ « 5 to A q = 5.9 comes from the phonon 
mode along F-N direction with q as mentioned above. At 
N symmetry point, the same mode contributes A q „ w 2.9 
to A q = 5.2. Smaller contributions with A q „ > 1 also 
come from q-points lying close to T— P direction and 
to- wards H point. The fact that the smaller regions in 
reciprocal space make substantial contributions to the 
electron-phonon coupling, a precise determination of A 
using a uniform q-grid becomes computationally diffi- 
cult. Since most of these contributions come from the 
low energy phonons, it is not expected to have a signifi- 
cant effect on the determination of T c . To illustrate the 
point, we also show in Fig. 6, the F(uj) and a 2 F(uj) of 
bcc yttrium calculated at the lattice constant a = 5.1987 
a.u. as described earlier. Not surprisingly, the difference 
in A, resulting from the two sets of calculations, is more 
pronounced in 15-25 meV region. 

The most important result of the present work is ob- 
tained by solving numerically the isotropic gap equation 
Til . T3| using the calculated Eliashberg function a 2 F(uj). 



The results of such a calculation for bcc yttrium as a 
function of [i* is shown in Fig. 7. We find that for fi* 
ranging from 0.1 to 0.25, the T c values range from 47.7 
to 36.8 K. A comparison of our calculated T c , includ- 
ing the ones described in Ref. [12|, with experiment is 
shown in Fig. 7. Thus, the present work clearly shows 
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Figure 6: The phonon density of states F(uj) (top panel) and 
the Eliashberg function o?F{ui) (bottom panel) of bcc (solid 
line) and hep (dashed line) yttrium. The dotted line in the 
two panels show F(lo) and a 2 F(u>) of bcc yttrium at a = 
5.1987 a.u., respectively. The two insets show u and A gj/ (for 
the lowest mode) along F-N direction. In the top panel, the 
phonon density of states of hep yttrium has been reduced by 
a factor of 4 for clarity. 
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Figure 7: The calculated T c (solid circles) as a function of /x* 
(left panel) and v/vo (right panel). Experimentally observed 
T c (open squares), taken from Ref. [3], are also shown in the 
right panel. The calculated T c for hep and dhep phases of 
yttrium are taken from Ref. [12]. 



the possibility of a T c of 30 K or above for compressed, 
bcc yttrium. 

In conclusion, we have shown that compressed, bcc 



yttrium has a large electron-phonon coupling with A = 
1.8, leading to a superconducting transition temperature 
T c k 30 K or above. We have also shown that such a large 
coupling is a direct consequence of lattice hardening and 
its subsequent coupling to Fermi electrons. 
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